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High-pr 7T° suppression in Au+Au collisions at -s/s^, = 200 GeV. 
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Neutral pions with px = 1 - 8 GeV/c have been measured for 9 centrality classes in 
Au+Au collisions at y/s^ = 200 GeV by the PHENIX experiment at RHIC. The vr° 
multiplicity in central reactions is significantly below the binary collision scaled yields 
from both peripheral Au+Au and pp reactions. The observed suppression sets in for the 
50-70% centrality class and increases with px and centrality. For the most central bin, 
the deficit amounts to a factor ~2.5 at p? ~ 2 GeV/c gradually increasing to a factor ~6 
at pt ~ 8 GeV/c. 

1. Introduction 

One of the primary goals of the PHENIX experiment at BNL RHIC is the measurement 
of hadron spectra out to large transverse momentum (p?) in high-energy heavy-ion (HI) 
collisions. High-pr particles result from the fragmentation of quarks and gluons produced 
in large Q 2 parton-parton scattering processes during the earliest stages of a HI collision 
and provide direct signatures of the partonic phase of the reaction. Interestingly, hard 
scattering multiplicities can be quantitatively compared (after scaling with the number 
of binary nucleon-nucleon, NN, inelastic collisions N co u) to: (i) baseline "vacuum" (pp) 
and "cold medium" (pA) data, and (ii) perturbative QCD predictions. In both cases, 
any departure (deficit or excess) from the "expected" results provides information on the 
strongly interacting hot and dense medium created in the reaction. One of the most 
intriguing results from the first RHIC run was the suppressed yield of moderately high- 
Pt tt° (pt ~ 1.5 - 4.0 GeV/c) in central Au+Au collisions at y/s NN = 130 GeV with 
respect to the scaled pp extrapolation [ [j]. Such a suppression is in qualitative agreement 
with theoretical predictions of parton energy loss effects ("jet quenching") in an opaque 
medium, the amount of energy loss being proportional to the reached gluon density [ |J. 



2. Experimental setup and data analysis 

During the 2001-2002 RHIC run, the full PHENIX Electromagnetic Calorimeter EMCal 
was instrumented providing a solid angle coverage of approximately Arj = 0.7 and A<f) = n. 
The integrated luminosity of y/s NN = 200 GeV Au+Au collisions was about 25 /ib _1 and 
the collected neutral pion statistics was a factor of ~ 10 2 larger than in Run-1, allowing 
the exploration of the 7r° spectra up to much larger values in pt and for finer centrality 
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classes. PHENIX also measured the 7r° spectrum in pp collisions at the same \/s [|3[]. The 
combination of full acceptance, high statistics, and the measurement of pp data in the 
same detector permits a very precise study of the high-p T n° suppression at 200 GeV. 

The data presented in this analysis correspond to pions from 30M AuAu "minimum 
bias" events with vertex position \z\ < 30 cm. The PHENIX EMCal consists of 4 lead- 
scintillator (PbSc) sectors (2592 towers per sector with 5.25x5.25x37.0 cm 3 size) in the 
west central arm plus 2 more sectors in the east one, and 2 lead-glass (PbGl) sectors (4608 
towers of 4.0x4.0x40.0 cm 3 ) in the east arm. The large radial distance of the calorimeter 
to the interaction region, of ~5 m, keeps the detector occupancy reasonably low (<15%) 
even in the highest multiplicity events of HI collisions at RHIC. Event centrality is deter- 
mined by correlating the charge and energy measured in two global detector systems: the 
Beam-Beam Counters and the Zero Degree Calorimeters. Neutral pions are reconstructed 
through an invariant mass analysis of 7 pairs detected in the 8 active EMCal sectors. The 
raw 7r° spectra are then corrected for geometric acceptance and (multiplicity-dependent) 
reconstruction and efficiency losses. The corrections were determined embedding 2M sim- 
ulated 7T° into real events. The final systematic errors in the fully corrected spectra are of 
the order of ±25%. Two independent analyses were carried out for the PbSc (shown here) 
and PbGl calorimeters yielding consistent results within the overall systematic error. 

3. Results 

Figure [1] shows the resulting 7r° pt distributions measured for 70-80% peripheral (left) 
and for 0-10% central (right) AuAu collisions, compared to the n° yield measured in pp 




Figure 1. Invariant tt° yields measured in peripheral (left) and in central (right) AuAu 
collisions (circles), compared to the N co u scaled pp 7r° yields (stars) [B|. The lines indicate 
the systematic errors in the scaled pp yields due to the uncertainties in N co u and luminosity. 



3 



collisions [ |3| scaled by the corresponding number of binary collisions given by a Glauber 
model calculation: N coU (70 - 80%) = 12.3 ± 4.0, and A^ co ^(0 - 10%) = 975 ± 94 respec- 
tively. Whereas ir° production in peripheral reactions is consistent with the point-like 
scaling expectation, the spectrum from the most central reactions is clearly suppressed 
(by a factor ~2.5 at ~ 2 GeV/c increasing up to a factor ~6 at p? ~ 8 GeV/c). A 
complementary way to depict such a suppression is using the nuclear modification factor 

Raa{pt) ~ (N coll )d 2 N;;/d v dp^ (1) 

which quantifies the deviation of the tt° yield in AA collisions with respect to the pp be- 
haviour (i.e. with respect to the absence of nuclear-medium effects) in terms of suppression 
or enhancement (Raa smaller or larger than unity, respectively). Figure 0a shows the pt 
dependence of Raa for n° emitted in 0-10% central AuAu reactions at two center-of-mass 
energies: 200 GeV {circles) and 130 GeV [ (stars). Three interesting conclusions can 
be extracted from this plot abovef^pr = 2 GeV/c: 



1. The Raa values at 200 GeV and 130 GeV are compatible (within their correspond- 
ing systematic errors) and clearly below unity: Raa(2 GeV/c) ~ 0.4 gradually 
decreasing down to Raa(8 GeV/c) ~ 0.16 for the highest y/s. 

2. The low Raa values measured at RHIC are clearly at variance with the high-p^ 
hadron production enhancement (Raa >1) observed at CERN-SPS energies (WA98 
data) [|,§: No "Cronin effect" [| is observed. 

3. The trend of Raa below unity in the range px = 2 — 8 GeV/c (corresponding to 
parton fractional momenta x = 2p T /yfs ~ 0.02 — 0.1) seems to be inconsistent 
with initial-state nuclear effects ("shadowing" of the Au parton distribution func- 
tion). Indeed, shadowing is known [ [/J to decrease for larger x (turning into "anti- 
shadowing" above x ~ 0.1), whereas Raa consistenly diminishes for larger pt values. 



The same quenching is evident in the nuclear modification factor constructed not from 
the ratio of AuAu to pp pions but from the ratio of central to peripheral AuAu tt° 
spectra (scaled by their corresponding centrality-dependent N coU ). Such Raa (see fig. 
5 of S. Mioduszewski contribution [ Q), -which have part of the systematic errors on 
the spectra cancelled out-, are also well below one in central collisions. Finally, fig. 
shows the evolution of the n° suppression at three fixed pr bins (pr= 2.2, 4.2 and 6.2 
GeV/c) as a function of the reaction centrality. The transition from the binary-scaling 
behaviour apparent in the peripheral region (Raa ~ 1, centrality class > 70%) to the 
highly suppressed central region is gradual and seems to occur over the 50-70% centrality 
bin (the systematic uncertainty in Raa, not shown in this plot, is of the order of ~20%). 

3 For pt < 2 GeV/c, Raa is below unity in all cases since the bulk of particle production is due to soft 
processes which scale with the number of participant nucleons (N part ) rather than with N co u. 
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Figure 2. Left. Nuclear modification factor for ir° emitted in top 10% central AuAu colli- 
sions at 200 GeV (circles), and 130 GeV (stars) (the bands show the syst. uncertainties). 
Right Raa at 200 GeV as a function of the reaction centrality for 3 fixed px bins. 



4. Conclusions 

Transverse momentum spectra of neutral pions have been measured at mid-rapidity 
by the PHENIX EMCal up to p T = 8 GeV/c for 9 different centrality classes from 30M 
Au+Au "minimum bias" events at y/s NN = 200 GeV. The spectral shape and invariant 
yield of peripheral reactions are consistent with that of pp collisions scaled by the number 
of inelastic NN collisions. Central yields, on the other hand, are significantly lower 
than peripheral AuAu and pp binary-scaled extrapolations, in agreement with the results 
found at y/s NN = 130 GeV. The observed quenching sets in over the 50-70% centrality 
class and increases with p? and centrality, being as high as a factor 6 at pt = 8 GeV/c in 
the top 10% central collisions. For these collisions, the nuclear modification factor Raa 
is systematically well below one and decreases gradually with p?, at variance with the 
common initial-state "phenomenology": (i) no "Cronin enhancement" is observed, and 
(ii) the growing suppression with px seems inconsistent with the decreasing dependence 
of "shadowing" on x. Since initial-state nuclear effects seem not determinant, final-state 
medium effects appear as a promising explanation for the observed high-pr pion deficit. 
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